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1 Introduction 

Short-rate models are of fundamental importance in the quantitative field of finance, 
as they provide a comprehensive mathematical framework for pricing interest rate or 
credit derivatives mm . The diversity of model structures and assumptions, enable 
us to choose the most appropriate approach when dealing with specific pricing issues. 
The basic Gaussian affine models, such as Vasicek lfl5l and Hull and White j6) gained 
interest among practitioners due to their analytical tractability and transparency, with 
closed-form pricing formulas available on hand. There is however, a trade-off be¬ 
tween such advantages and implausible model forecasts, which allow negative in¬ 
terest rates. Some others, such as the Cox, Ingresoll and Ross model J4j, despite 
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having the property of positive rates, often provide unrealistic outcomes by imposing 
a lower positive bound for the par swap rate m . Non-negativity of interest rates does 
not seem so important, or may even be undesirable in today’s low-interest rate envi¬ 
ronment. However, in the context of default intensity modeling, negative hazard rates 
are generally not feasible due to lack of consistency with arbitrage-free assumptions. 

A model developed by F. Black and P. Karasinski (BK) in 1991, also known as 
the ’’exponential Vasicek model” a, overcomes the problem of negative rates. It 
postulates log-normality of short rates, motivated by the fact that the market stan¬ 
dard Black formulae for caps and swaptions are based on log-normal distributions 
of relevant rates. Moreover, it possesses rather good fit-to-data properties, especially 
concerning the swaption volatility surface. Unfortunately, in this model, exact ana¬ 
lytical formulae for swaptions, or even for zero-coupon bond prices, do not exist. 
This lack of analytical tractability requires the use of computationally intensive and 
time-absorbing numerical methods (PDE or Monte Carlo). This virtually precludes 
efficient model calibration and seriously narrows areas of potential model applica¬ 
tions. 

In response to challenges related to implementation of the BK model, there have 
been several attempts to obtain reliable analytical approximations of zero-coupon 
bond or swaption prices. In particular, Tourrucoo et al. fT4l proposed approximate 
formulas for zero-coupon bonds in a one-factor model, derived in the limit of small 
volatility by applying the regular asymptotic expansion of a transformed bond PDE 
fll . Antonov and Spector |2| went further and came up with a generalised multi-factor 
BK model. By performing a regular asymptotic expansion of the PDE, they provide 
approximations for both zero-coupon bond and European swaption prices. Never¬ 
theless, both approaches consider small volatility cases, which are not plausible in 
many financial applications, in particular in credit markets, where lognormal volatil¬ 
ities of default intensity as high as 100% are frequently observed ED- A different 
approach towards deriving approximations has been adopted by Stehlikova IfTZj . who 
developed small time expansions for one factor models. Prices of zero-coupon bonds 
were represented by means of a Taylor series expansion with coefficients represented 
in a closed form, obtained via specific recurrent relations. Another approximation 
concept, which originated from chemical physics under the name of ’’the exponent 
expansion”, was introduced to Finance by Capriotti 0 and applied to the calculation 
of transition probabilities and Arrow-Debreu prices, for several diffusion processes. 
This approach appeared to provide very accurate approximations and was further pur¬ 
sued by Stehlikova, Capriotti 02 in the context of the BK model. On the basis of 
the exponent expansion, the authors proposed representing the price of a zero-coupon 
bond in the form of a power series in time, that can easily be computed by means of a 
recursion involving only simple one-dimensional integrals. For larger time horizons, 
the exponent expansion can be combined with a fast numerical convolution to obtain 
more accurate results. However, these approximation formulas are quite complex and 
there is generally no guidance on the selection of the proper truncation order of the 
exponent expansion, to obtain a sufficient degree of approximation accuracy. 

In this paper, we propose a totally different approach towards approximating zero- 
coupon bond and swaption prices. The concept is based on a novel technique, apply¬ 
ing the Karhunen-Loeve representation |8] of the Ornstein-Uhlenbeck process and 
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one other related process, that appear in the context of a BK model. Initially, we pro¬ 
vide semi-analytic approximations of zero-coupon bond prices. Afterwards, we de¬ 
rive analogous approximations for swaptions using a conceptually similar, yet more 
elaborate approach. The formulae are easy to implement, computationally fast and 
provide very accurate approximations for the vast majority of parameter settings. As 
zero-coupon bonds and swaptions (caps/floors as a special case) are basic instruments 
used for calibration of the model, our results can be used directly for that purpose, 
substantially improving the speed of calibration. 

The structure of this paper is organized as follows. In the second section we intro¬ 
duce the mathematical framework established for deriving price approximations in a 
BK model. The third section provides two kinds of approximations of a zero-coupon 
bond, with all related derivations included. The fourth section is composed in a simi¬ 
lar manner and deals with swaption pricing. The final fifth section contains numerical 
results of approximations for bonds and swaptions for several selected sets of model 
parameters, as well as comparisons with alternative approximations for bond prices 
based on papers lfT4l . lfm . For the readers convenience, most technical and purely 
mathematical considerations are gathered in the Appendix. 

2 Mathematical preliminaries 

We denote by (r t ) t >o the process of the short rate and by /, := ln(r,) it’s natural log¬ 
arithm. We assume that the process (/;) f >o follows dynamics postulated in a Black- 
Karasinski model 0i-e. 


dl t = ( a(t) — bit) dt + odW t , 


( 2 . 1 ) 


where a,b are positive constants, a(t) is some deterministic function of time and 
(W f ) ; >o is a Wiener process under the spot measure. Let us introduce some notations 
and recall a few facts. 

Notation: For « > 0, v > 0 we set 



U 



(2.3) 


X, (u) := —(/„+,-e b 'l u -A(u,u + t)). 

o 


(2.4) 


In the specific case u = 0 we denote X t := Xf°\ 


Proposition 1 For any u,t > 0 


r u +t = r u/ exp(t7X, (,,) ) 


(2.5) 
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where >_ is an Ornstein-Uhlenbeckprocess satisfying 

X 0 (h) = 0, dX t (u] = -bX t {u) dt + dW r iu) , (2.6) 

where Wf" 1 W u + t — W u is a Wiener process. 

Proof Using (12. Il l and elementary calculus we check that (x/“^V >0 is an Ornstein- 
Uhlenbeck process satisfying (12.6I >. By virtue of introduced notations we have 

r u+t = exp (A(u,u + t) + e~ bt \n(r u ) + aX^ u) \ = r uf ex p (ox/"*). (2.7) 

□ 


Remark 1 Processes (W^) ;>0 and (x/ 1 ^) >( , are independent of JF„. 

Remark 2 Quantity r u v can be interpreted as the value of the short term rate r u+v in 
the absence of volatility. It is also the dominant of the distribution of r u+v under ■ C P U . 

Remark 3 Ornstein-Uhlenbeck process (X f ) r >o is a centred Gaussian process with a 
covariance function 0 

K(sf) := 4 e “ &|t “ S| - i e “ 6(l+S) (2-8) 

2b 2b 

and a variance function 

V(s):=X( S ,s) = ^(l- e “ 2fo ) (2.9) 

For further purpose we introduce yet another process 
Notation: Let us define an Ornstein-Uhlenbeck bridge process T as 

X s :=X s -^22x t , J€[0,r], (2.10) 


Remark 4 Ornstein-Uhlenbeck bridge (X t ) t ^ Q r j is a centred Gaussian process satis¬ 
fying X(): Xj 0, with a covariance function 


K(sf) := K(s,t) 


K(s,T)K(t,T) 

V(T) 


and a variance function 


V(s) := K(s,s)=V(s) 


e -2b(T-s) 

l-e- 2bT 


( 2 . 11 ) 


( 2 . 12 ) 


Moreover for any s G [0.7’] a random vector (Xj,X s ) has a joint Gaussian distri¬ 
bution and by definition E(X 7 ’X J ) = 0, hence X s is independent of Xj. 
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Now we can state main theorems, on which our approximations are based on. 
Here we present their assertions only. For detailed proofs please refer to the Ap¬ 
pendix. 

Theorem 1 (Karhunen-Loeve Theorem, Loeve ||8l ) Let (X t ) te ^ ab j be a centred 
stochastic process with a covariance function K. ThenX, admits the expansion (called 
the Karhunen-Loeve representation or Karhunen-Loeve expansion) 


X t = Y, \[Kfn(t)Z n a.s., 

n =0 


(2.13) 


where 



(2.14) 


are orthogonal random variables such thatEZ n = 0 andEZ^ = 1. Moreover {f n {t )} n >0 
form an orthonormal basis in Lr([a,b]) consisting of the eigenfunctions (correspond¬ 
ing to non-zero eigenvalues) of a Fredholm operator SZ , i.e. 


ZZ fn = fnfn, n>0. (2.15) 

The series (12. 1 3b converges a.s. and uniformly for t £ [a,b] in the norm ||F(/)|| := 

1/9 

(ET 2 (f)) . In particular, if(X t ) te \ a ^ is a centred Gaussian process, then Z n ’s given 

by (12.14b are independentN (0,1) random variables. 


The following two theorems give the explicit form of Karhunen-Loeve expansion 
for the Ornstein-Uhlenbeck process and Ornstein-Uhlenbeck bridge. 


Theorem 2 Let (X t ) te ^ibe an Ornstein-Uhlenbeck process in the interval [0,t] 
satisfying the equation (12.6b . Then its Karhunen-Loeve expansion is of the form 


X,= Y, (2.16) 

72=0 

where Z n are independent, identically distributed normal N(0 ,1) random variables 
and are given by 




where 0]„(z) is the unique solution of the equation 

cocot(ft)T) = —b 

in the interval f (n + l)f j 5 n £ NU{0}. 


(2.17) 


(2.18) 
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Theorem 3 Let {X t ) ; g o, 7 ’] be an Ornstein-Uhlenbeck bridge process in the interval 
[0,7’] satisfying (12.10b . Then its Karhunen-Loeve expansion is of the form 

£= £ \Jln(T)f nJ (t)Z n , (2.19) 

77=1 

where Z„ are independent, identically distributed normal N(0, 1) random variables 
and f n j,kn{T) are given by 

/»jW= = <2 - 0) 

For the purpose of further elaborations we also introduce a few important nota¬ 
tions and recall some useful theorems and lemmas. 


Notation: Let H n denote the probabilistic Hermite polynomial of degree n and let 
h„i,..., h nu denote its zeros (in ascending order). 

Notation: For a function / : R —> R let 2„ (/) denote its Lagrange interpolation 
polynomial with nodes li n i.e. 


lUf)} (z) := I f{h n ,k) / _ , , t • 

k=l (Z nn,k)Wji=k\finJi bn.j) 


Notation: For a function / : R —)• R let Q„ (/) denote the Gauss-Hermite quadrature 
of degree n, rescaled to the probabilistic convention, with abscisas h n j. and weights 


w n ,k, i-e. 




&«(/) : = Y, W n,kf(h n ,k), 

(2.21) 


k= 1 


where 

2 n ~ 1 n ! 




(2.22) 


Remark 5 It is well known that for functions 
cally holds 



lim Q„(/) = 

ft_y_]_oo 



dz, 


i^dx'j asymptoti- 
(2.23) 


which becomes an equation for polynomials of degree up to 2n — 1. 


Notation: Define the operator § : R[X] —t R[X] and the functional h : R[X] —> R such 
that for polynomial W (z) = LaLo w k? k 

77 — 1 

[W] (z) := J>*z*. $(W) :=«_!, 

A—0 


(2.24) 
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where coefficients u ,t are defined recursively as follows 

u n = 0, u n+ 1 = 0, (2.25) 

Mi = (k + 2)u k+2 + w k+ i, k = n- 1,...,-1. (2.26) 


Remark 6 Note that §. I) are properly defined without the implicit assumption deg(W ) = 
n. Indeed, for any m > n setting 0 = w m = ... = w n+ 1 and W* (z) = T!k=o w kZ k we ob¬ 
tain 

500 = 60*0, f>(w*) = t)(w). (2.27) 

Lemma 1 For any polynomial W (z) 

- Jw(z)e-^dz=[§>{W)](z)e-^-\){W)0{z)+C. (2.28) 

For the proof refer to Appendix. 

Lemma 2 For any continuous, bounded function f : R — > R 

lim £„(/)=/, (2.29) 

n —>+<*> 

where convergence holds in the L 2 norm. 

For the proof refer to the Appendix. 


3 Approximation for zero-coupon bonds 

In this section we present two semi-analytical formulae for approximate pricing of 
zero-coupon bonds in a BK model. The derivation of the first formula is based on 
general pricing fundamentals and application of Theorem|2] The second formula adds 
some additional approximations and simplifications, making it more usable without 
much loss of accuracy. 


Notation: Let T, z > 0,n > 0. We define the following functions 


Fn,i(t,zo, ■ ■ ■ ,Zn) := exp <7 £ \/%(f)fk,x(t)Zk 


k =0 


Gn,x(t) ■■= exp (v(t) - 52 4(t)A,t( 0 2 ) j » 

T 

/„ jT (zo,...,z„) := J r T ,tG n . T {t)F r , iT (f ,z 0 ,... ,z«) dt. 


(3.1) 


(3.2) 


(3.3) 
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Approximation 1 The sequence 

(-“> +' 

1 

B 


„(r,r+T):-(^_) n+i |."|e x p( ~° + 2 + ~") 


provides approximation 


B(T,T + t)= lim fl„(7\r + T). 

n —f-°° 


rfzo • • • 

(3.4) 

(3.5) 


Proof Considering the process (x/ r M as defined in Proposition|T|and it’s Karhunen- 

V /fe[0,r] 


Loeve expansion from Theorem|2] we have 

rr+t = Ttj ex P ( (JX t {T) ) = r* n T+t R n ^(t), 


where 


r n,T+t — r T fG n ^(f)F n ^(t ,Zo,... ,Z„), 

£h,t( 0 = G,7,T( f ) ex P ( E a Vh(T:)fk,x{t)Zk ) • 

\fc=n+1 / 


(3.6) 

(3.7) 

(3.8) 


Notice that rj t is a deterministic, continuous function of t, r* T+t is -independent 

and ^ 

Er* r+( =r T , t exp , r€[0,T],n>0, (3.9) 

ER BiT (0 = 1, t € [0, x),n > 0, (3.10) 

E< T (f)=exp((J 2 E 4 WA,t( 0 2 J, <e[0,T],«>0. (3.11) 

V k=n +1 / 

Fix f S [0, t]. Recalling formulas from Theorem|2] we notice that 


4t 2 


sup A, t (0 < k k (T) < 2 2 , 

f e[o,T] T (2k+l)-n 2 


therefore from (13. 1 11 1 

ER 2 x (t) < exp 
Hence by Jensen’s inequality 

E|R„,t(0 - 1| < \J E (R n ,r( l 

Finally, since 


8t <t_ +? 1 

k 2 Jth ( 2 * +1 ) 2 


= :M„. 


(3.12) 


(3.13) 


= ^er2,(0- 

1 < 

/ 5 

\ 

-y< r -K* 

V 0 



(3.15) 
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then applying the inequality |exp(— x) — exp(—y)| < \x—y\ for x,y > 0, switching the 
order of integration and using (13.9b . (13. 14b we obtain 


B{T, T + t) - B n (T, T + t) | < E ^ J (r T+t - r* nT+t ) dt 

1 \ T 

j\ r T+i ~ r *nj+t | dt J = j ^r* lT+t E\R n , x {t)-\\dt 

o Jo 

<T sup ?T.t exp f ^—V(T) \ 

refO,Tl V ^ / 


(3.16) 


and since M„ —> 0, we get our assertion. 


□ 


Remark 7 Although the formula from Approximation[T]may seem quite complicated, 
it is applicable and for small n can easily be calculated numerically. In particular, the 
external integrals can be evaluated with high accuracy by the use of the (probabilistic) 
Gauss-Hermite quadrature, whereas the internal integral I n T can be calculated using 
the Romberg method or Legendre quadrature (provided that a(t) satisfies sufficient 
smoothness conditions). 

Remark 8 ApproximationQ]presents the convergence property of the sequence 
{B„(T,T + f)}„ '"()■ For practical applications however, one shall use approximation 
B„(r, T + T) for a specific integer n. Fortunately, the first approximation (n = 0) al¬ 
ready appears to give very accurate results. This can be understood by taking into 
account that the first eigenvalue in the Karhunen-Loeve expansion accounts for the 
dominant part of the overall variance of the Ornstein-Uhlenbeck process. Based on 
this observation, we introduce a new series of zero-coupon bond price approxima¬ 
tions. 


Approximation 2 The sequence 

+ t) := &„(exp(— 7 0i t)), n> 1 (3.17) 

provides approximation 

B 0 (T,T+t)= lim B* n (T,T + x). (3.18) 

n —>+°° 

Proof Note that Iq t (z) is a strictly positive, continuous function of z £ R. Therefore 
exp(—7 q T (-)) G L 2 f R, -^=e~T x ~dxj and by (12.23b we get our assertion. 

Important property of such approximation is its monotonic dependence on rj. 

Proposition 2 For any n > 1, B*(T,T + t) is a strictly decreasing function of rj. 
Moreover 


lim B*(T,T + T) = 1, lim B*(T,T + t)=0. 

rj —>0 rr-r+oo 


(3.19) 
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Proof By definition 

B*(T,T + t) = Q„(exp(—/ 0 , T )) = £ w n,k^p(~kA h n,k))- (3.20) 

k= 1 

Since weights w n ^ are positive, then it is sufficient to prove that lo.rihn.k) is a strictly 
increasing function of r-/ for each k. Recalling the formulae 

T 

Io,i( h n,k) = J r T jGo,T{t)Fo, x {t,h nt k)dt, 1 <k<n, (3.21) 

o 

we note that Go,r,Fo. T are strictly positive functions independent of rj, whereas fy, 
is a strictly increasing function of rj. Consequently, /o. T (/i„i) is a strictly increasing 
function of ry. 

Moreover, since Go **(■)> Fo.kx ('>hn, j) are bounded on [T,T + t] and convergences 


lim rj t = 0, lim rjt — +°° 

rj —>0 ’ rj —>+oo 


are uniform in t £ [T,T + t], we obtain 

lim exp (-/ 0 jvt(/u-j)) = 1 , lim ex P ( ~ h),Ni{h,j)) = 0- 

rj —>0 rj —>+od 


Since by definition 


B*(T,T +k8) 


Lw„,;exp( 
J =1 


and that w n j sum up to 1, we get our assertion. 


(3.22) 


(3.23) 


(3.24) 

□ 


4 Approximation for swaptions 

Let us consider the swaption with expiry T > 0, strike S > 0 and underlying swap 
of tenor T = /Vi), where n> 1 is the number of fixed-leg payments and 5 > 0 is the 
length of payment period. Consider the parameter ft) related to the swaption type, 
equal ft) = 1 for a payer swaption and ft) = -1 for a receiver swaption. 


Notation: Denote the stochastic discount factor and its conditional expectation in 
respect to Xj by 


P(t) :=exp j r s dsj , j3(x) :=B(f}(T)\X T =x). (4.1) 

Denote the values of swap annuity, underlying swap rate and two auxiliary quan¬ 
tities as 

N 

A(t,T,N) := 5 ^ B(t,T +kS), 

k= t 


(4.2) 
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B(t,T)-B(t,T+NS) 

{ ’ ’ h A(t,T,N) 

C(T,N,S) :=B(T,T +N8)+SA(T,T 1 N) 1 

P(T,N,S;X t ):=P(X t )(i-C{T,N,S\X t )^, 


(4.3) 

(4.4) 

(4.5) 


where the last argument Xj expresses (implicit) dependence of the variable on Xj 
and can be omitted when not relevant. 

Proposition 3 The swaption price Swpt can be expressed by the formula 

Swpt = ^(g)^-{(o>wC(t,n,s)}P{T,N ,S',Xt)^ . (4.6) 

Proof Pricing the swaption under a spot measure as the expectation of its intrinsic 
value at expiry and using our notations (14.3I) - (14.4I) . we can express the theoretical 
price of a payer/receiver swaption as 

Swpt = E(/3 (T)A(T, T,N ) (a(r(T, T,N) - s) )' + ) = E (j3 (T) (to (1 - C(T,N, S )) ) + 

(4.7) 

By the tower property, this expectation can be calculated taking conditional expec¬ 
tation with respect to Xj first. Since C(T,N,S ) is ^-measurable this leads to our 
assertion. 

□ 

To obtain the approximation for swaption price we will apply approximations to 
ft (Xt) and C(T,N,S) in (14.6I >. Let us start with approximating P{Xt). 


Notation: Let m > 1. We define the following functions 

( m r - ^ is(f r r\ \ 

\Jh{T)f k , r {t)zk + o V ( T J xJ : 

G m ,z{t) ■= exp L^( r )A,r(t) 2 )^ , 

1 

Imj(x,Zl,Zm) '■= J rT,tG m ^(t)F m ^(t,X,Zl, ■ ■ ■ ,z m ))dt. 


(4.8) 

(4.9) 

(4.10) 


Approximation 3 The sequence 

Pm{x) ■= J ■ ■ ■ J eX P , • • • ,Z m ) - + ) dzi ... dz„ 


(4.11) 


P(x)= lim P, n (x). 

m—>-+°° 


provides approximation 


(4.12) 
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Proof The proof is analogous to that of Approximation!]] modified so that all random 
variables are replaced with their expectations conditional on Xj. All functions and 
constants are replaced with their corresponding ”hat” versions, dependent on x. Such 
analogy is based on equation (12.10b and the fact the Ornstein-Uhlenbeck bridge and 
all Zfc are independent of Xj. 

□ 

Approximation 4 The sequence 

j3*(x) := Q,„(exp (-I m y)(x,-)), m> 1 (4.13) 

approximates j3o(x) 

A)W= lim /3*(x). (4.14) 

m—>+oo 

Proof The proof is a direct analogue of that presented in Approximation^ 

□ 

Notation: For m.n > 1 let us denote 

P min (T,N,S)=P,UT,N,S-X T ) := fa(X T ) (l - C* n (T,N,S',X T )), (4.15) 

N 

C*{T,N,S-Xt) := B*(T,T+N8) + S8 £ B* n {T,T+k8), (4.16) 

k= 1 

where the last argument Xj expresses (implicit) dependence on AY, which is clear 
from the fact that each B* n (T, T + kS) is a function of ly and rj = ro.rexp^AY) (ro.r 
is a constant). 

Proposition 4 For any m,n > 1 there exist xo = xoijn.n) such that I\ n n (T,N,S\ ■) is 
negative on (—°°,xo) and positive on (xo,+°°). 

Proof Note that rj is an increasing function of xj and 

lim rj = 0, lim rj = +°°. (4.17) 

Therefore using Proposition [2] we observe that C*(T,N,S;x) is a strictly decreasing 
function of Xj and 

^lim^C*(r,A,5;x) = 0, tim^C*(T,N,S;x) = 1+N8S >\. (4.18) 

Hence there exists finite xq such that 

C*(T,N,S-x)> 1, x<x 0 and C*(T,N,S'x)< 1, x>x 0 , (4.19) 

which concerning I5*(Xt) > 0, is equivalent to our assertion. 

□ 

Proposition 5 For any m.n > 1 the price Swpt of the swaption can be approximated 
with the formula 

Swpt « Swpt* n n E(^(ot[ co x T >cox 0 }Pm,n(T,N,S',Xr)j . 


(4.20) 
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Proof Approximations!]]-EJdirectly imply that 

N 

P{T,N,S\X t ) « Po(X T )(l -B 0 (T , T +N8) - S8 £ B 0 (T, T + kS)) 

k= l (4.21) 

«^(X r )(l-C;(T,lV,5;X r ))= J P m] „(r,lV ) 5;X r ). 

In addition, from Proposition [4] we deduce that the condition (oP{T,n,S\Xj) > 0 is 
equivalent to coXt > ftUo, hence finally 

Swpt « E^l {fflXr > (nA . o} P ; , lj „(7’,A^,5;X r )), (4.22) 

which is our assertion 

□ 

Now we are ready to provide a tractable formula for approximate swaption pric¬ 
ing. To this end we simply replace Swpt,* n in the last approximation with its specific 
Lagrange interpolating polynomials. 

Approximations Form,n,k> 1 denote 

fm,n ■ fm,n (z) = P m ,n{T,N,S ; y/V(T)z), (4.23) 

= Wn,n). (4.24) 

Take any k,l > 1 such that 

ftn.n (h k .l ) < 0 < fn,n{fk,l+ 1) (4.25) 

andzt E (hk,hhic t i+i\ such that 

fm,n,k(4) = 0- (4-26) 

77zen f/ze sequence 

SwpC^k := ^ { [$(£**)] (4)exp (4.27) 

provides approximation 

Swpt* = lim Swpt* k . (4.28) 

’ *->+<» 

In particular, the swaption price Swpt can be approximated as 

Swpt » (4.29) 

Proof First of all note that required k,l,zt exist. Indeed 

— h k .\ = h k ,k ->• +°°, £ -> +°°. (4.30) 

Hence by Proposition^ 

fm,n(h k , l) < 0 < fm,n{h k ,k) (4.31) 


for almost all k and 


”4 0 < fm,n(h k ,l+ 1) 


(4.32) 
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for a specific /. Refering to the definition of /* ., this also means 

/m,n(M <°</m,«(%,/+!)> (4.33) 

hence f^ n k has a root z* k in {hu,h u +\]- 

By Proposition 0 and Lemma[T|we have respectively 

SwpC„ = E (cot {B)Z > mo} f(Z )), (4.34) 

ft)oo 

S™pt* m , n , k =~^ j = E (o)1 {£oZ > q , 4} ^„ ] ,(Z)) , (4.35) 

Z k 

where Z ~ A(0,1), zo : = ^ • Then denoting 

£k ■— E ^l{coZ>a)zo} — /m.«(Z)^ , (4.36) 

4 := E (ail {oZ > a4} (/ ra ,„(Z) (Z))) , (4.37) 

we have 

Swpt* nn - Swpt* n n k = £ k + 8 k . (4.38) 

We will show that £ k . 8 k —> 0. First note that 

|P( 0 )Z> o)zo)-P(toZ> 0JZ?)| = |®(zo)-«P(zJ)| < -^=|zo-z*|. (4.39) 

Additionally, from elaborations in Proposition [4] 

\Pm,n(T,N,S;X T )\ <max{l,A55'} < 1 + N8S, (4.40) 

hence f m < 1 +N8S. Combining those together 

\e k \<-j=(l+N8S)\zo-z* k \. (4.41) 

By definition of xo from Proposition [H we note that zo £ (h k .i,h k i + ]]. Taking in ac¬ 
count that z* k also lays in this interval, we have 

\e k \ < -jL(l +N8S)(h ktl+l -h ktl ) (4.42) 

and noting that h k / + \ — h k j —> 0 as k —> +°° (see Theorem 6.1.2 na) we get e k —t> 0. 
Now, let us observe that 

8 2 <E\f{Z)-f^ nk (Z)\ < ||/-/;^||, (4.43) 

where ||-|| denotes the norm in L 2 ^R, — L=e~? x2 dx^J. As f, kn is a continuous and 
bounded function, thus by Lemma[2]it holds ||/* k — f m . n \\ —> 0 as k —> +°°, hence 
8 k —>■ 0. Finally, since £ k + 8 k —»• 0, then taking into account Proposition0and (14.381) . 
we finally get our assertion. 
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□ 

Remark 9 Despite some complexity in the approximation formulas obtained, they 
are easily computable in practice. Namely, proper / can be found by evaluating 
P* n (T,N,S; y/V ( T)hj ) at zeros of Hermite polynomial (j = 1,... ,k), which can be 
easily done numerically. Those same values are also used for calculating the coeffi¬ 
cients of polynomial f* k , specifically by solving a set of linear equations. After¬ 
wards, coefficients of the polynomial § {f * t n .) and the value I)(/^ „ ,) are directly 
computable by using recursive formulas given by their definition. Finally, the root 
z* k of f* n n k can be effectively found with any standard numerical procedure (such 
as Newton’s method or the false position method). The only numerically extensive 
element here is the evaluation of P* in , which requires k(Nn + m) numerical integra¬ 
tions ( k arguments hj, N maturities of bonds approximations, each consisting of n 
components of quadrature, plus m nodes of quadrature for approximation of j3). 

Remark 10 For practical applications, one should take some specific values. 

In typical situations we recommend using k = m = n = 5, which proved to be ac¬ 
curate enough in our numerical tests and require a moderate number of numerical 
integrations. 

Remark 11 Note that for exact quantities it holds 

E(j3(X r )(l—C(7>,S))) =l-C(0,n,S). (4.44) 

However, it does not remain true if we substitute p (AY) and C(T.n.S) with their 
approximations applied in Approximation [5] Consequently, our approximations do 
not obey the put-call parity exactly. 


5 Numerical results 

In order to test the accuracy of the approximations presented in previous sections, 
several prices of zero coupon bonds have been computed. For simplicity a(-) was as¬ 
sumed to be a constant function of time of the form a(t ) = Fln(r avg ), where r avg = 3%. 
For better clarity and comparability, results are presented in the form of yields-to- 
maturity (with a continuous compounding convention), not actual prices of bonds. 
In order to examine the dependence of results on different parameters, yields-to- 
maturity were calculated for the following: 

• maturities: 1, 2, 5, 10 and 20 years 

• values of ro: 1%, 3% and 6% 

• values of b: 0.02 and 0.1 

• values of cr: 25% and 50% 

In Table 1 we present results of such calculations, obtained from Approximation 
13 benchmarked to the exact results obtained via Monte-Carlo simulations. One can 
see, the errors of approximations are very small, in the order of at most a few basis 
points for every set of parameters examined. 
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Similar numerical tests were performed for swaption prices. Several prices of 
payer and receiver swaptions have been computed, for various sets of parameters of 
the model and the swaption payoff. We used the same sets of model parameters as for 
bonds. For each of them, we examined the accuracy of approximations for various 
swaptions, using parameters as follows: 

• swaption expiries: 1, 2, 5 and 10 years, 

• underlying tenors: 1, 2, 5 and 10 years, 

• swaption moneyness: 10%, 80%, 90%, 100%, 110%, 125%, 150%. 

(moneyness being a quotient of the swaption Strike and ATM strike as the forward 
swap rate). Due to the dimensionality of the parameter space, we present results in 2 
layers: 

• dependence on moneyness for a given tenor (Table 2), 

• dependence on underlying tenor for ATM strike (Table 3). 

In the case of strikes other than ATM, only out-of-the-money swaptions were con¬ 
cerned, i.e. the price of payer or receiver swaption was calculated depending on 
whether its strike was above or below the forward swap rate (such restriction can 
be imposed without loss of generality, because of the put-call parity). 

As in the case of bonds, prices calculated using Approximation0were compared 
to ’’exact” prices (calculated on the lattice). In addition to such checks, we examined 
the scale of put-call disparity resulting from approximation, as mentioned above in 
Remark Qj] To this end, in Table 4 we compared prices of ATM-payer and ATM- 
receiver swaptions (which in principle should be equal) obtained from Approxima¬ 
tion [3 To allow for easier comparison of results for different swaptions, and to stay 
compatible with market conventions, all swaption prices were translated into their 
implied volatilities. Hence figures in tables represent differences between implied 
volatilities corresponding to compared swaption prices (approximate vs exact in the 
case of Tables 2 and 3 or payer vs receiver in Table 4). 

As you can see, with only a few exceptions, errors of approximations range from 
-50 bp to +50 bp, where the vast majority are less than 10 bp in terms of absolute 
value, which is far below a typical bid-offer spread. Not surprisingly, the biggest 
errors are observed in the case of long expiry/tenor and/or high volatility. Similar ob¬ 
servations address the put-call disparity, which appears negligible except in cases of 
the longest expiries and/or tenors, which reflect weaker efficiency of approximations 
when applied to longer time horizons. 

Finally, we compared our approximations for bonds with those obtained via meth¬ 
ods proposed in . We used Table 1 from mu herein, which contains approxi¬ 

mations of zero-coupon bond prices obtained with both those methods as well as from 
Monte Carlo simulations, as benchmark values. Calculations have been conducted for 
specific sets of parameters, namely: 

• maturities: 0.1, 0.5, 1, 2 and 3 years, 

• r 0 = 6%,b = ln(0.04), cr = 85%. 

Table 5 includes these results complemented with prices obtained from our Ap¬ 
proximation [2] However, in order to keep our convention, we converted bond prices 
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to their yields and present the results in the form of a difference vs benchmark (MC). 
The results in Table 5 reveal moderately good performance of Approximation [2] in 
comparison to other approaches. Most importantly, it maintains a rather stable er¬ 
ror rate while increasing bond maturity, whereas results from Ref [16] exhibit very 
different behaviour, with errors increasing strongly with maturity. 



18 


A. Daniluk, R. Muchorski 


Table 1 Yield-to-maturities obtained by Monte Carlo simulations (MC) and Approximation |2](A2). 


Model Parameters 

Maturity 


Yield 


?0 

b 

a 

MC 

A2 

Error 




i 

1.071 % 

1.071 % 

0.000% 




2 

1.142% 

1.142% 

0.000% 

1% 

0.1 

25% 

5 

1.350% 

1.350% 

0.000% 




10 

1.663 % 

1.663 % 

0.000% 




20 

0.832 % 

0.832 % 

0.000% 




1 

3.043 % 

3.043 % 

0.000% 




2 

3.080 % 

3.080 % 

0.000% 

3% 

0.1 

25% 

5 

3.159 % 

3.159 % 

0.000% 




10 

3.221 % 

3.222 % 

0.001% 




20 

1.610% 

1.611 % 

0.001% 




1 

5.885 % 

5.885 % 

0.000% 




2 

5.769 % 

5.769 % 

0.000% 

6% 

0.1 

25% 

5 

5.435 % 

5.436 % 

0.001% 




10 

4.971 % 

4.975 % 

0.004% 




20 

2.485 % 

2.487 % 

0.002% 




1 

1.027 % 

1.027 % 

0.000% 




2 

1.053 % 

1.053 % 

0.000% 

1% 

0.02 

25% 

5 

1.134% 

1.134% 

0.000% 




10 

1.264% 

1.264% 

0.000% 




20 

0.632 % 

0.632 % 

0.000% 




1 

3.046 % 

3.046 % 

0.000% 




2 

3.089 % 

3.089 % 

0.000% 

3% 

0.02 

25% 

5 

3.203 % 

3.203 % 

0.000% 




10 

3.331 % 

3.333 % 

0.001% 




20 

1.666 % 

1.666 % 

0.001% 




1 

6.048 % 

6.048 % 

0.000% 




2 

6.086 % 

6.086 % 

0.000% 

6% 

0.02 

25% 

5 

6.145 % 

6.146% 

0.001% 




10 

6.075 % 

6.081 % 

0.006% 




20 

3.038 % 

3.041 % 

0.003% 




1 

1.120% 

1.120% 

0.000% 




2 

1.243 % 

1.243 % 

0.000% 

1% 

0.1 

50% 

5 

1.607 % 

1.607 % 

0.000% 




10 

2.104% 

2.107% 

0.003% 




20 

1.052% 

1.053 % 

0.001% 




1 

3.178 % 

3.178 % 

0.000% 




2 

3.336 % 

3.336 % 

0.000% 

3% 

0.1 

50% 

5 

3.668 % 

3.670 % 

0.002% 




10 

3.872 % 

3.882 % 

0.009% 




20 

1.936% 

1.941 % 

0.005% 




1 

6.137 % 

6.137 % 

0.000% 




2 

6.215 % 

6.216 % 

0.001% 

6% 

0.1 

50% 

5 

6.174% 

6.181 % 

0.006% 




10 

5.747 % 

5.774 % 

0.026% 




20 

2.874 % 

2.887 % 

0.013% 
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Table 2 Differences between implied volatilities corresponding to swaptions prices calculated using Ap¬ 
proximation [5]and pricing on lattice, calculated for various moneyness levels. 


Model Parameters 


Swaotion 



Implied volatility error vs 

moneyness (%ATMF) 










RECEIVER 



PAYER 


ro 

b 

<7 

Mat. 

Ten. 

Fwd IRS 

70% 

80% 

90% 

100% 

100% 

110% 

125% 

150% 




1 

1 

1.22% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

1.0% 

0.10 

25% 

2 

2 

1.43% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

5 

5 

1.99% 

0.00% 

0.00% 

0.01% 

0.01% 

0.01% 

0.01% 

0.00% 

0.00% 




10 

10 

2.60% 

-0.01% 

0.00% 

0.01% 

0.01% 

0.02% 

0.02% 

0.01% 

0.00% 




1 

1 

3.17% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

3.0% 

0.10 

25% 

2 

2 

3.25% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

5 

5 

3.34% 

-0.01% 

0.00% 

0.01% 

0.01% 

0.02% 

0.02% 

0.01% 

-0.01% 




10 

10 

3.31% 

-0.01% 

0.00% 

0.00% 

0.00% 

0.03% 

0.03% 

0.02% 

0.01% 




1 

1 

5.82% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

6.0% 

0.10 

25% 

2 

2 

5.46% 

-0.01% 

0.00% 

0.00% 

0.00% 

0.01% 

0.01% 

0.00% 

0.00% 

5 

5 

4.62% 

-0.01% 

0.00% 

0.01% 

0.01% 

0.03% 

0.02% 

0.01% 

0.00% 




10 

10 

3.86% 

-0.01% 

-0.01% 

-0.01% 

-0.01% 

0.03% 

0.03% 

0.03% 

0.02% 




1 

1 

1.09% 

0.00% 

0.00% 

-0.01% 

-0.01% 

0.00% 

0.00% 

0.00% 

0.00% 

1.0% 

0.02 

25% 

2 

2 

1.17% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

5 

5 

1.40% 

0.01% 

0.02% 

0.02% 

0.01% 

0.02% 

0.01% 

0.00% 

-0.01% 




10 

10 

1.72% 

0.11% 

0.12% 

0.11% 

0.10% 

0.12% 

0.09% 

0.05% 

-0.01% 




1 

1 

3.18% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

3.0% 

0.02 

25% 

2 

2 

3.30% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

5 

5 

3.52% 

0.03% 

0.05% 

0.05% 

0.04% 

0.06% 

0.05% 

0.02% 

-0.02% 




10 

10 

3.51% 

0.01% 

0.02% 

0.02% 

0.02% 

0.07% 

0.06% 

0.06% 

0.05% 




1 

1 

6.32% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

0.00% 

6.0% 

0.02 

25% 

2 

2 

6.38% 

0.00% 

0.00% 

0.01% 

0.01% 

0.01% 

0.01% 

0.00% 

-0.01% 

5 

5 

6.20% 

0.01% 

0.03% 

0.03% 

0.03% 

0.06% 

0.05% 

0.04% 

0.01% 




10 

10 

5.34% 

-0.07% 

-0.13% 

-0.17% 

-0.18% 

-0.13% 

-0.12% 

-0.10% 

-0.02% 




1 

1 

1.38% 

-0.02% 

-0.03% 

-0.04% 

-0.03% 

-0.02% 

-0.01% 

0.00% 

0.02% 

1.0% 

0.10 

50% 

2 

2 

1.75% 

-0.02% 

-0.02% 

-0.02% 

-0.01% 

-0.01% 

-0.01% 

0.00% 

0.01% 

5 

5 

2.63% 

0.23% 

0.23% 

0.21% 

0.17% 

0.23% 

0.18% 

0.10% 

-0.02% 




10 

10 

3.26% 

0.04% 

0.08% 

0.09% 

0.07% 

0.27% 

0.25% 

0.21% 

0.14% 




1 

1 

3.55% 

0.00% 

-0.01% 

-0.01% 

-0.01% 

0.00% 

0.00% 

0.00% 

0.01% 

3.0% 

0.10 

50% 

2 

2 

3.90% 

0.07% 

0.08% 

0.07% 

0.06% 

0.08% 

0.06% 

0.02% 

-0.03% 

5 

5 

4.16% 

0.23% 

0.26% 

0.25% 

0.22% 

0.32% 

0.27% 

0.20% 

0.07% 




10 

10 

3.93% 

-0.06% 

-0.06% 

-0.07% 

-0.08% 

0.16% 

0.15% 

0.16% 

0.18% 




1 

1 

6.49% 

0.01% 

0.02% 

0.02% 

0.01% 

0.03% 

0.02% 

0.01% 

-0.01% 

6.0% 

0.10 

50% 

2 

2 

6.43% 

0.12% 

0.15% 

0.14% 

0.12% 

0.15% 

0.12% 

0.06% 

-0.03% 

5 

5 

5.49% 

0.10% 

0.13% 

0.13% 

0.12% 

0.27% 

0.25% 

0.21% 

0.15% 




10 

10 

4.38% 

-0.10% 

-0.17% 

-0.22% 

-0.24% 

0.01% 

0.01% 

0.04% 

0.13% 
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Table 3 Differences between implied volatilities corresponding to swaptions prices calculated using Ap¬ 
proximation [5] and pricing on lattice, calculated for various swaptions expiries and underlying tenors, for 
ATM strike. 


Model Parameters 

Tenor 

ATMF Payer Volatility Error 





Expiry 





ro 

b 

O 

1Y 

2Y 

5Y 

10Y 




1 

0.00% 

0.00% 

0.00% 

-0.01% 

1 % 

0.1 

25% 

2 

0.00% 

0.00% 

0.00% - 

-0.01% 

5 

0.00% 

0.00% 

0.01% 

0.00% 




10 

0.01% 

0.02% 

0.01% 

0.01% 




1 

0.00% 

0.00% 

-0.01% - 

-0.03% 

3 % 

0.1 

25% 

2 

0.00% 

0.00% 

0.00% - 

-0.02% 

5 

0.01% 

0.01% 

0.01% 

0.00% 




10 

0.02% 

0.02% 

0.01% 

0.00% 




1 

0.00% 

0.00% 

-0.02% 

-0.06% 


0.1 

25% 

2 

0.00% 

0.00% 

-0.01% - 

-0.04% 


5 

0.02% 

0.02% 

0.01% - 

-0.01% 




10 

0.02% 

0.02% 

0.01% - 

-0.01% 




1 

-0.01% 

0.00% 

-0.01% - 

-0.02% 

1 % 

0.02 

25% 

2 

-0.01% 

0.00% 

-0.01% 

0.00% 

5 

-0.01% 

0.00% 

0.01% 

0.02% 




10 

0.03% 

0.05% 

0.09% 

0.10% 




1 

0.00% 

0.00% 

-0.01% 

-0.04% 

3 % 

0.02 

25% 

2 

0.00% 

0.00% 

0.00% - 

-0.02% 

5 

0.02% 

0.03% 

0.04% 

0.02% 




10 

0.14% 

0.15% 

0.12% 

0.02% 




1 

0.00% 

0.00% 

-0.02% 

-0.08% 

6% 

0.02 

25% 

2 

0.00% 

0.01% 

0.00% - 

-0.05% 

5 

0.05% 

0.06% 

0.03% - 

-0.04% 




10 

0.11% 

0.08% 

-0.05% - 

-0.18% 




1 

-0.03% 

-0.02% 

-0.01% - 

-0.13% 

1% 

0.1 

50% 

2 

-0.05% 

-0.01% 

0.01% - 

-0.07% 

5 

0.06% 

0.14% 

0.17% 

0.04% 




10 

0.46% 

0.49% 

0.32% 

0.07% 




1 

-0.01% 

0.01% 

-0.02% 

-0.20% 

3% 

0.1 

50% 

2 

0.02% 

0.06% 

0.04% - 

-0.10% 

5 

0.31% 

0.35% 

0.22% 

0.01% 




10 

0.58% 

0.49% 

0.15% - 

-0.08% 




1 

0.01% 

0.02% 

-0.05% - 

-0.33% 

6% 

0.1 

50% 

2 

0.09% 

0.12% 

0.03% - 

-0.18% 

5 

0.45% 

0.40% 

0.12% - 

-0.10% 




10 

0.46% 

0.29% 

-0.10% - 

-0.24% 
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Table 4 Call-put disparity for corresponding ATM payer and receiver swaptions evaluated using Approxi- 
mation[5] expressed as the difference of corresponding implied volatilities of payer and receiver swaptions. 


Model Parameters 

Tenor 

Payer-Receiver ATMF Volatility 




Expiry 





ro 

b 

O 

1Y 

2Y 

5Y 

10Y 




i 

0 .00% 

0 .00% 

0 .01% 

0 .02% 

1 % 

0.1 

25% 

2 

0 .00% 

0 .00% 

0 .01% 

0 .02% 

5 

0 .00% 

0 .00% 

0 .00% 

0 .02% 




10 

0 .00% 

0 .00% 

0 .01% 

0 .02% 




1 

0 .00% 

0 .00% 

0 .02% 

0.08% 

3% 

0.1 

25% 

2 

0 .00% 

0 .00% 

0 .02% 

0.06% 

5 

0 .00% 

0 .00% 

0 .01% 

0.04% 




10 

0 .00% 

0 .00% 

0 .01% 

0.03% 




1 

0 .00% 

0 .01% 

0.05% 

0.18% 


0.1 

25% 

2 

0 .00% 

0 .01% 

0.04% 

0 .12% 


5 

0 .00% 

0 .00% 

0 .02% 

0.07% 




10 

0 .00% 

0 .00% 

0 .01% 

0.04% 




1 

0 .01% 

-0.01% 

0 .01% 

0.04% 

1 % 

0.02 

25% 

2 

0 .00% 

0 .00% 

0 .01% 

0.03% 

5 

0 .00% 

0 .00% 

0 .00% 

0 .02% 




10 

0 .00% 

0 .00% 

0 .00% 

0 .02% 




1 

0 .00% 

0 .00% 

0.03% 

0 .11% 

3% 

0.02 

25% 

2 

0 .00% 

0 .00% 

0 .02% 

0.08% 

5 

0 .00% 

0 .00% 

0 .02% 

0.06% 




10 

-0.01% 

-0.01% 

0 .01% 

0.05% 




1 

0 .00% 

0 .01% 

0.06% 

0 .22% 

6 % 

0.02 

25% 

2 

0 .00% 

0 .01% 

0.04% 

0.16% 

5 

0 .00% 

0 .00% 

0.03% 

0 .10% 




10 

-0.01% 

-0.01% 

0 .01% 

0.06% 




1 

0 .01% 

0 .01% 

0.06% 

0.35% 

1 % 

0.1 

50% 

2 

0 .00% 

0 .00% 

0.06% 

0.28% 

5 

-0.01% 

0 .00% 

0.06% 

0.25% 




10 

-0.15% 

-0.10% 

0 .01% 

0.19% 




1 

0 .00% 

0 .02% 

0.15% 

0 .66% 


0.1 

50% 

2 

0 .01% 

0 .02% 

0.13% 

0.52% 


5 

-0.03% 

-0.01% 

0 .10% 

0.38% 




10 

-0.16% 

-0.09% 

0 .02% 

0.23% 




1 

0 .01% 

0.05% 

0.29% 

1 .10% 


0.1 

50% 

2 

0 .01% 

0.04% 

0.23% 

0.81% 


5 

-0.04% 

0 .00% 

0.15% 

0.49% 




10 

-0.14% 

-0.08% 

0 .02% 

0.25% 
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Table 5 Yield-to-maturities of bonds obtained from various approximations, expressed and an error vs 
Monte Carlo simulations (MC). The respective approximations are derived from small volatility expansion 
(Ref 1141 ). the exponent expansion truncated to the first EE(1), second EE(2) and third term EE(3) and our 
Approximation [2] 


Ref 03) vs MC 

EE(1) vs MC 

EE(2) vs MC 

EE(3) vs MC 

A2 vs MC 

-0.10% 

0 .00% 

0 .00% 

0 .00% 

-0.02% 

-0.23% 

0 .02% 

0 .02% 

0 .02% 

-0.01% 

-0.46% 

0 .01% 

0 .00% 

0 .00% 

-0.07% 

-0.90% 

0.06% 

0.03% 

0 .00% 

-0.13% 

-1.24% 

0.17% 

0 .10% 

0 .00% 

-0.08% 


Appendix 


Proof of Lemma [l] 

Proof At first note that operator § and functional l) are linear. Namely, let a be any 
real number and W\(z) = ffl=[) w \.kZ k ; W^fz) = L<” =0 w 2 ,kZ k polynomials. By Remark 
5, without loss of generality we assume they have the same degree. Denote 


[S(Wi)](z) = I «i,*z*, [S(w 2 )](z) = I «2 , k z\ 

k =0 k =0 

[$(«Wt)] (z) = £ 8 kf, [S(Wi + W 2 )\ (z) = £ h k z k , 

k—0 k =0 

then by definition 

U] j; = = U 2 = g/i = = = ^/z+1 “ 0^ 


(5.1) 

(5.2) 


(5.3) 


in specific 

hn = U \ n U2.m hn+l — tfljH-fl “t” t/2,n+li = tXU |/j, gl,«-|-l ~ ^ w l,n+l ■ (5.4) 


In addition, for k: — \ <k<n — 1 


gfc+2 = «Wl,fc+2> A 1^+2 = Ml,Ar+2 + M2,ifc+2) ^k+2 = Ml,i+1 + M2,t+1 > (5-5) 

thus 

gt = (k + 2)gA+2 + awi i A + t = a((* + 2)gi,*+2 + w i,*+i) = awi,A, (5.6) 

hk = (^ + 2)/tA+2 + M ; l,/t+l + Vl^2,it+1 ^ 

= ((£ + 2)mia+ 2 + Wl.tr+l) + ((k + 2)u2,k+2 + w 2,k+l) = Ml,* + M2,A, 

hence by the induction step it is clear that linearity properties are satisfied. 

Now define operator : R[x] —> R[x] and functional b* : R[x] — > R such that for any 
polynomial W (z) 

- J W(z)e~? z2 dz= [?>*(W)}(z)e-? z2 ~l)*(W)<P(z)+C. (5.8) 
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Let us note that 5v and b* are properly defined and linear, thus the lemma postulates 

s*=s, r=h, (5.9) 

hence by linearity of §*, h, h* it is sufficient to prove 

r(z") = 6(z"), W) = S(z") (5.10) 

for any non-negative integer n. To this end observe that the recursion formula for 
coefficients of § (z n ) provides 


(n - 1)!! 

u k = {(n — k) mod 2}———, 0 <k<n, 

m_ i = {(n + 1) mod 2} (n — 1)!!, 
whereas from integration by parts 

= z"- l e~ lZ +(n— l)z"~ 3 e~ lZ —(« — !)(« — 3) Jz n ~ A e~^ z dz = ... 


(5.11) 

(5.12) 


/'L ^ 1 ( n - l) m 

= j Y, {(« — k) mod 2}- 
\k =0 


it!! 


1 -2 


{(« + 1) mot/ 2}(n — l)!!<J>(z) +C, 

(5.13) 


hence coefficients of match respective values m* and h*(z") = m 


□ 


Theorem 4 (Nevai (9j) Let / /;c a continuous function defined on the real line and 
L n (f.x) the Lagrange interpolation polynomial interpolating f at zeros of the noil- 
probabilistic Hermite polynomial of degree n. Assume that f satisfies 

lim f(x)(l + \x\)e~^ = 0. (5.14) 


Then 


lim 

«-> + »O 



—oo 


L n (f,x) \e~& 


p 

dx = 0 


holds for every p > 1. 


Proof of Lemma [2] 


(5.15) 


Proof Consider the function 

/: fix) = / ( V2x) . (5.16) 

Since zeros of the probabilistic Hermite polynomial are s/2 times the corresponding 
zeros of the non-probabilistic Hermite polynomial (of the same degree), then setting 
Z = \/2z we have the following identity 

Ln(f,Z) = [£„(/)] (Z), 


(5.17) 
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where £„ is specified in Theorem[4] By changing the variable we have 

+°° ^ +°° 

J (/(z)-[ fi n(/)](z)) e~^'dz=y/2 J (/(z) -L n (/,£)) e~ zr dl (5.18) 

Notice that since / is a continuous and globally bounded function, then so is /. It’s 
clear that any bounded function satisfies condition <15.14b . hence setting p = 2 Theo- 
rem [4] postulates 

+o° 

e ~ ;2 dz = 0, (5.19) 

thus (15.181) implies our assertion. 

Lemma 3 Let K\,K 2 : [0, T} 2 -> K . be C l functions and also let f : [0, T] —> R and 
K : [0, T} 2 —> R be continuous functions, such that 

K(s,t) = K\(s,t), s<t, (5.20) 

K(s,t) = K 2 (s,t), s>t. (5.21) 

Denote 

• function K* f : (7f*/)(s) := fif K(s,t)f(t)dt 

• functions h' s ,h" s as the first and second-order partial derivative of h(s,t) with 
respect to variable s for any twice differentiable function h(s,t) 

• single argument function ~jli(s,s) as the first-order derivative of h(s,s) any dif¬ 
ferentiable single argument function li(s,s) 

Then we have the following: 

a) K* f is differentiable on [0,7'] and 

(K * f)'(s) = (K 2 (s, s) - K\ (s,s))f{s) + (K' s * /) (s) . (5.22) 

b) Moreover if K\ and K 2 belong to C 2 ((0,7) 2 ) and 

Ki{t,t)=K 2 (t,t) for all te[0,7], (5.23) 

then K * / is twice differentiable and 

(K*f)"(s)= ^^ 2 (s, S )-^i( S ,^/( S ) + (^"*/)( S ). (5.24) 

Proof We have 

5 T 

(- K*f)(s ) = j K 2 (s,t)f(t)dt + j Ki{s,t)f{t)dt (5.25) 

0 s 
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Hence by differentiating integrals we get 

S S 

J ' K 2 (s,t)f(t)dt = K 2 {s,s)f{s) + J K' s (s,t)f(t)dt (5.26) 

o o 

and 

T T 

^jK l (s,t)f{t)dt = -K 1 (s,s)f{s) + j K' s (s,t)f(t)dt. (5.27) 

S S 

Summing up those equations we get the first part of our assertion. Assuming that 
K\ (s,s) = K 2 (s,s), equation (15.22b obtains the form 

(K*f)\s) = (K' s *f)(s), (5.28) 

which can be rewritten 

5 T s T 

jK 2 (s,t)f(t)dt + ^~ JK\(s,t)f{t)dt = J'K' s (s,t)f(t)dt + JK' s (s,t)f(t)dt. 

0 i 0 i 

(5.29) 

Differentiating it once more and replacing in (15.22b function K' s (s,t) with K" (s.t) we 
finally get (15.24b . 

□ 

Corollary 1 Let K and K be covariance functions from (12.8b and (12.11b respectively 
and f be a continuous function on [0,7’]. Then we have the following 

a) Function K * f is differentiable on [0,7] and 

{K*f)'{s) = -b[K*f){s). (5.30) 

b) Functions K* f,K * f are twice differentiable on (0,7) and 

(K*f)"(s)= ^K 2 (s,s)-^K l (s,s) S jf(s) + (K'; s *f)(s), (5.31) 

(K*f)"(s)= (J-K 2 (s,s )-/(*) + (*> /)(*)• (5-32) 

Proof Notice that functions K. K are jointly continuous functions such that K can be 
represented by functions K\ for s <t and K 2 for s > t, where K\,K 2 are of the form 

K X (sf) = 1 (e~ b ^ - e~ b ^) , (5.33) 

K 2 (s,t) = ± (e- b(s ~ t} - e - fe(,+i) ) ( 5 . 34 ) 
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and function K can be represented by functions K\ for s < t and K 2 for s > t, where 
K \, K 2 are of the form 


Ki{s,t) =Ki{s,t) 


2bK(s,T)K(t,T) 

1 — e~ 2bT ’ 


(5.35) 


K 2 (s,t) = K 2 {s,t) 


2bK(s,T)K(t,T) 
1 — e~ 2bT 


(5.36) 


Moreover 

K l {t, t ) = ^(\- e - lb ') =K 2 (t,t ) (5.37) 

and 

*i(M) = Yb 0 ~ e ~ 2bt ) - T_l-2?r (5.38) 

thus assumptions in Lemma [3] are met for functions K\ ,K 2 ,K and K\. K 2 . K. Partial 
derivatives of functions K\,K 2 satisfy 


Y s k YY = ~bKj(s,t), j = 1,2, (5.39) 

hence from Lemmaj^we have the first part of the assertion. 

Now, computing first and second derivatives provides 



s) 


d_ 

Ts 


K\ (s,s) = (-b-b)- 


J_ -6(s-s) 

2b 


- 1 , 


(5.40) 


d » <9 ~ d d 

—K 2 (s,s) - -^Ki{s,s) = -^K 2 (s,s)- -^Ki{s,s) = -1 

and 

d 2 

-^Kjisj) = b~Kj(s,t), j = 1,2, 

d 2 „ 

j^Kjisj) =b 2 Kj(s,t), j = 1,2. 

Therefore applying Lemma [3] we obtain 

(K*f)"(s) = -f(s)+b 2 (K*f)(s) 

and 

(- K*f)"( S ) = -f( S )+b 2 (K*f )( S ), 

which is the second part of our assertion. 


(5.41) 

(5.42) 

(5.43) 

(5.44) 

(5.45) 

□ 


Proof of Theorem |3 
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Proof Due to Theorem |T| we know that the Ornstein-Uhlenbeck process (X,) fg [ 0 T | 
has an expansion of the form (12. 1 6b . where Z„ are independent /V(0.1) random vari¬ 
ables. Hence it is sufficient to obtain closed formulas for A„(t) and /„. x . Let K be 
the covariance function of the process (X, ) /e - 0 r as in (12.8I >. We proceed with finding 
eigenfunctions of the operator associated with kernel K. Recalling (12.151) we know 
that functions /„. x satisfy the equation 

Sf fn,T = (5.46) 

for some respective A„(t) >0 and by Lemma [3] we can alternatively write 


K* fn,x — ^n{y)fn,x- (5.47) 

Differentiating both sides twice and making use of Corollary[T|we get 

(5.48) 

Therefore 


-fnM +b 2 (K*f n , T ) (s) = A„(t)/''». 
~fn,x{s) +b 2 ?i n (l)fnj(s) = A„(t)/" t (s) 


or equivalently 


(A„(t)Z? 2 - 1 )/ n , T (s) = A„(t)/" (s), 


(5.49) 


(5.50) 


hence the solution is a linear combination of functions exp(±zft)„f), where ft),, is such 
that 


2 1-A„(t )b 2 

co,: = 


A„(t) 


(5.51) 


Moreover function f n _ T satisfies two boundary conditions. Firstly, we notice that 
K(0,t) = 0 for t £ [0, t], therefore 


i 

0) = J K(0,t)f n , r (t)dt =0. 


(5.52) 


Secondly, applying Corollary Q] we observe that 

An(T)/' ]T (T) = (K*f n>x )\ T) = -b(K s *f n , x ){T) = -R,(t)/„ iT (t). (5.53) 

Since ft),, can be real or pure imaginary, depending on the sign of the right-hand side 
of (15. 51k let us consider those two cases in detail: 


1° (On £ iR. (O n 4 0 

In this case ft),, = ia> n for some ft),, £ R\ {0} and - is a linear combination of func¬ 
tions exp(±ft)„f). Hence considering condition (15.52b . we deduce that f n T is of the 
form 

fn,r(t) =c„(exp(—ftj„f)-exp(<a„f)), (5.54) 

where c„ £ R\ {0} (since /„ T is a real, not identically equal zero function). Thus 
applying condition (15.53b to function /„ T , we obtain from (15.54b the equation 

ft)„(exp(—ft>„T) +exp(ft)„T)) = fr(exp(—ft>„T) — exp(<B„T)), 


(5.55) 
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which cannot be satisfied since the left and right-hand side of the equation have op¬ 
posite signs (as b > 0). Consequently, this case provides no solutions. 

2° (On £ R 

In this case f n T is a linear combination of functions exp(±/fi)„f), so condition (15.521) 
implies that 

fn.i(t) = c' n (exp(i(O n t) — exp(— i(0„t)) = 2ic' n sin(ov) =c„sin(o„f), (5.56) 

where co n 4 0 and c n = — \ic' n for some c' n £ R\ {0} (since f n r is a real, not identically 
equal zero function). Hence relying on the form of /„ T from (15.561 ). condition (15.531) 
provides 

G)„cos(o)„t) = — bsm((o„T), (5.57) 

which after elementary calculations can be rewritten as 

B,{(O n x) =-bx. (5.58) 

Setting c, (x) :=xcotx. Equation (15.571 ) implies that sin(e)„T) 4 0, thus function ^ is 
well-defined and continuous on each interval («7T, (n 4-1 )tt), n £ NU {0}. Moreover 
it is strictly decreasing on each such interval. Indeed 

, x sin(2x) — 2x 

| (t) = cotx-=---< 0. (5.59) 

sin x 2sin“x 

Therefore for each 77 € NU |0} function C„ c \ / , is invertible. One can 

also easily see that (n + 1)7Z^ = (—°°,0), hence equation (15.581 ) has for 

each n £ N U {0} exactly one solution of the form 

co„ := -^~\-bx). (5.60) 

X 

The value of A„(t) corresponding to each (0„, implied from (15.5 1 b . is given by 

W := I- (5 ' 61 > 

b z + (Of[ 

thus this case provides solutions to (15.501) . that satisfy conditions (15.521) . ( 15.531 ) and 
have the general form 

fn,r(t) =c n sin((0„t), A„ = -,-1— 2 n = 1,2,... (5.62) 

b A + (On 

Now, in order to prove that functions f n T and corresponding A„(t) are respectively 
the eigenfunctions and eigenvalues satisfying (15.471) . for each n £ NU {0} we define 
function g n such that 


g„(t) = (K*f n )(t)-K{i:)fn^{t) 


( 5 . 63 ) 
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and prove that it is identically equal zero. First of all, from the properties of f n x and 
A,,(t), we can easily deduce that 

g'hG) = {K*f n ,r)"(,s) - A„(T)/" T (i) 

o / \ o o (5.64) 

= (-fn,x{s) + b 2 {K*f ntX ) (s)) - (-/„,t(s) +^A„(t)/„, t (s)) = 

and because function g n is real, then the solution to (15.64b is of the form 

g„(s) = d\ exp(/?f) +a? 2 exp(— bt), c/i^GK.. (5.65) 

Moreover we have 


gni 0) = (K * /„ jT ) (0) — A„(t)/„ iT (0) = 0 (5.66) 

and 

SnM = t (t) 

(5.6/) 

= -Z?(fi:*/ niT )(T)+H l (T)/„ ;T (T) = -bg n {x). 

Equation (15.661 1 implies that c /2 = — c/| , hence from (15.671 ) we have 

c/i(exp(/?T) +exp(—fer)) = — c/i(exp(Z?T) — exp(— bx)), (5.68) 


which holds if and only if d\ =0, because for non-zero d\ the left and right-hand side 
of the equation have opposite signs (since b > 0). Thus d\ = di = 0 and functions g„ 
are identically equal zero, as desired. 

Finally, to obtain f UT normalized (as required in the Karhunen-Loeve expansion) we 
find the appropriate c n by the direct calculation: 


sin(2a) = 


1 = ll/n|| 2 = J c 2 sin 2 (co„t) dt 
0 

Taking into account that 
and that by (15.571 ) 

we have 

which leads to the equation 


= c„ 


/ T sin(2fi)„T) 
4 (O n 


2 tan a 
1 + tan 2 a 


1 =ci - + - 


r) = 

dif* 

I 

II 


2 ©„/? 


b 2 + co 2 ’ 

1 

b \ 

2 

b 2 + co 2 J 


hence 


Gn — 


1 \ -1/2 

-(t + bX n (T)) 


T + bXnpy 


(5.69) 

(5.70) 

(5.71) 

(5.72) 

(5.73) 

(5.74) 


Substituting into (15.621) values c„ and A,„(t) from (15.741 ) and (15.611 ) respectively, we 
finally get the assertion. 
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□ 

Proof of Theorem U 

Proof Due to TheoremQ]we know that the Ornstein-Uhlenbeck bridge process (X; ) /<E [o. 7 ’] 
has an expansion of the form ( 12.19b , where Z n are independent /V(0.1) random vari¬ 
ables. Hence it is sufficient to obtain closed formulas for A n {T) and f n j. Let K be the 
covariance function of the process (Xf)f<=[o,r] as i n <12.1 lb . We proceed with finding 
eigenfunctions of the operator Sf associated with kernel K. Recalling (12.15b we know 
that functions f n j satisfy the equation 


^ fn.T = UT)fn,T 


(5.75) 


for some respective A n (T) > 0 and by Lemma[3]we can alternatively write 



(k*f„. t)(s) = XfT)f n j{s). 

(5.76) 

Differentiating 

both sides twice and making use of Corollary Q] we get 



- fn,T {S) +b 2 (K* f nJ ) ( S ) = (i) ■ 

(5.77) 

Therefore 

-/»,r(i) +b 2 X n (T)f ntT (s) = A „(T)?„' iT (s) 

(5.78) 

or equivalently 

( l n {T)b 2 -\)f nJ {s) = l n {T)f'; j {s ), 

(5.79) 

hence the solution is a linear combination of functions exp(±/ft),,f), where ft),, 
that 

&2 1 -\,(T)b 2 

" A„(r) 

is such 

(5.80) 


Moreover function f rLT satisfies two boundary conditions. Firstly, we notice that 
K(0,t) = 0 for t £ [0, T], therefore 


T 

i(7)/„,r(0) = J K{Q,t)f n ;r(t)dt = 0. (5.81) 

o 

Secondly, we observe that K(T,t) = 0 for t £ [0, T], thus 

T 

ln(T)fn, T (T) = J k{T,t)f n ;r(t)dt = 0. (5.82) 

0 

Since a>„ can be real or pure imaginary, depending on the sign of the right-hand side 
of (15.80b . let us consider those two cases in detail 


1° d)„ £ iR, d>„ + 0 
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In this case d) n = ico* for some co* £ R\ {0} and f nx is a linear combination of func¬ 
tions exp(±co*l). Hence considering condition (15.8 1 b . we deduce that f n j is of the 
form 

= c n (exp{-co*t) -exp(( 0 *t)), (5.83) 

where c„ £ R\{0} (since f n j~ is a real, not identically equal zero function). Thus 
applying condition (15.82b to f n .r, we obtain from (15.83b the equation 


(exp(-0)„*r)-exp«r)) =0, 


(5.84) 


which cannot be satisfied unless c„ = 0 or co* = 0. Consequently, this case provides 
no solutions. 

2° co f ] £ R 

In this case f n j is a linear combination of functions exp(±/o>„f), so condition (15.811) 
implies that 


f„j(t) = cj, (exp(/a)„f) — exp(—idv)) = 2ic' n sm{& n t) = c„sin(a>„f), (5.85) 


where c„ = — \ic' n for some c' n £ R\ {0} (since f n j is a real function). Hence, relying 
on the form of f n j from (15. 81b , condition (15.82b provides 


c„sin(d) ; ,7’) = 0 

has for each n £ N exactly one solution of the form 

nTTt 



(5.86) 


(5.87) 


The value of X n (T) corresponding to each ft)„, implied from (15.51b . is given by 



(5.88) 


thus this case provides solutions to (15.79b that satisfy conditions (15.81b . (15.82b and 
have the general form 



Now, in order to prove that functions f n j and corresponding X n {T) are respectively 


the eigenfunctions and eigenvalues satisfying (15.76b . for each n £ N we define func¬ 
tion g„ such that 


gn(t) = (K * fn tT ) (t) ~Xi(T)fn,T 


(5.90) 


and prove that it is identically equal zero. From the properties of f n j and we 

can easily deduce that 


§:(*) = (K*fn, T )"(s)-X n (T)% T (s) 

= (-fn,T(s) + b 2 (K*f n j)(s))-(-f n j(s)+b 2 X n (T)f n j(s))=b 2 g n (s) 


(5.91) 
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and because function g n is a real function, then the solution to (15.911) is of the form 
g n {s) = d\ exp(Z?f) +c/ 2 exp(— bt), di,d2&^. (5.92) 

In addition we have 

&(0) = (K * fnj) (0) — %n{T)f n j(Q) = 0 (5.93) 

and 

gn(T) = (K*f n j){T)-%n(T)f n j(T)=0. (5.94) 

Equation ( 15.93b implies that d2 = —d\, hence from ( 15.94b we have 

d i (exp(fcr) - exp(-for)) = 0, (5.95) 


which holds if and only if d\ = 0 (since b > 0). Thus d\ = d 2 = 0 and functions g n 
are identically equal zero, as desired. 

Finally, to obtain f n T normalized (as required in the Karhunen-Loeve expansion) we 
find the appropriate c n by the direct calculation 


1 = 



(5.96) 


hence 



(5.97) 


Substituting into ( 15.89b values c„ and X n (T) from (15.97b and (15.88b respectively, we 
finally get the assertion. 


□ 
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